Multicriticality and entanglement in the one-dimensional quantum compass model 
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We study the one-dimensional (ID) quantum compass model with two independent parameters 
by means of an exact mapping to the quantum Ising model. This allows us to uncover hidden 
features of the quantum phase transition in the ordinary one-parameter quantum compass model, 
showing that it occurs at a multicritical point where a line of first-order transitions intersects a 
line of second-order symmetry-breaking transitions of Ising type. We calculate the concurrence and 
the block entanglement entropy in the four ground state phases, and find that these entanglement 
measures accurately signal the second-order, but not the first-order, transitions. 
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I. INTRODUCTION 

The orbital degeneracy in d-shell transition metal ox- 
ides underpins many of the fascinating phenomena found 
in these materials. Well-known examples include ferro- 
electricity, colossal magnetoresistance, and charge order- 
ing [l[ . An early attempt to model the directional nature 
of the orbital states in the case of a two-fold degeneracy 
was made by Kugel and Khomskii 0], who introduced a 
simplified model — the quantum compass model (QCM) 
— where the orbital degrees of freedom are represented 
by (pseudo)spin-l/2 operators and coupled anisotropi- 
cally in such a way as to mimic the competition between 
orbital orderings in different directions. For example, 
on a two-dimensional (2D) square lattice, the coupling 
between spin components on neighboring lattice sites is 
tied to the direction of the corresponding bonds, with 
Ising-like interactions J^crf CTf_^_- and JyO-JaJj^- along the 
X- and y-axes of the lattice. Much interest has focused on 
possible "order-from-disorder" phenomena in this model. 
While the competition between spin interactions produce 
a massively degenerate ground state at the classical level, 
the degeneracy gets lifted by thermal or quantum effects, 
favoring a directional ordering of spin fluctuations at low 
temperatures 0, 01 • Added interest in the 2D model has 
been spurred by the proposal that it describes the physics 
of a collectively generated "protected qubit", with a su- 
perconducting Josephson junction array as a possible re- 
alization Q. Also, the model has been shown to be dual 
to the Xu-Moore plaquette model oip-\-ip superconduct- 
ing arrays and more recently, to Kitaev's toric code 
model in a transverse magnetic field 

The degeneracies in the energy spectrum of the quan- 
tum compass and related orbital models make numerical 
simulations very demanding. In the absence of exact so- 
lutions in two or higher dimensions, progress in mapping 
out the phase diagram of the model is coming only slowly 
and piecewise. There is strong evidence for the existence 
of a symmetry-broken ground state in 2D Q, however, 
the character of the quantum phase transition (QPT) 
into this state has been a controversial issue. A contin- 
uous second-order transition is favored in one st udy [8 1 , 
with others supporting a first-order transition [lol . uA \ . 



Recently Brzezicki et al. [l^ considered a one- 
dimensional (ID) version of the quantum compass model, 
with the aim of analytically exploring how the interplay 
between quantum fiuctuations and competing spin inter- 
actions produces a QPT. While the ID model may not 
directly bear upon the physics of its higher-dimensional 
relative, it still serves as an interesting counterpoint. By 
a clever construction, where pairs of spin operators on 
neighboring sites are altcrnatingly mapped onto terms 
in a quantum Ising model (QIM), Brzezicki et al. [l^ 
obtained an exact solution for the ground state energy, 
revealing that the ID model exhibits a first-order transi- 
tion between two disordered phases with opposite signs of 
certain local spin correlators. Intriguingiy, this first-order 
transition was found to be accompanied by a diverging 
correlation length for spin correlations on one sublattice. 

In this paper we exploit the method of Brzezicki et al. 
[l^ to study an extended version of the ID QCM, ob- 
tained by introducing one more tunable parameter. The 
expanded parameter space allows us to uncover the true 
character of the QPT identified in Ref. [H. As follows 
from our analysis, this transition in fact occurs at a mul- 
ticritcal point where a line of first-order transitions meets 
with a line of second-order transitions. The second-order 
critical points are Ising-like and separate locally ordered 
ground states from disordered states (with ordered quasi- 
local correlations). In contrast, when going through one 
of the first-order transitions, the only noticeable change 
in the character of the ground state is that pair correla- 
tions of spins on every second bond flip sign. In an effort 
to probe the anatomy of the various ground states we 
have analytically calculated the concurrence for pairs of 
spins, and also the block entanglement entropy. The con- 
currence signals the second-order QPTs, and the block 
entropy shows that they are Ising-like. None of these en- 
tanglement measures is affected by the first-order transi- 
tions. 

The rest of the paper is organized as follows: In Sec.HIl 
we introduce the extended ID QCM, and apply the 
method of Brzezicki et al. [l^ to obtain its spectrum 
and spin correlators. In Sec. IIIII we exploit our exact re- 
sults to map out the zero-temperature phase diagram of 
the model, and identify the character of its QPTs. In 
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Sec. IIVI we then perform an entanglement diagnostic of 
the various phases and phase transitions. Sec. |Vl finally, 
contains some concluding remarks. 

II. THE ONE-DIMENSIONAL QUANTUM 
COMPASS MODEL 

Consider the Hamiltonian 

N' 

H = ^[Jl(^2i-1^2i + J'2.'^2i-1^2i + Lia2i(J2i+i ], (1) 
i=l 

with periodic boundary conditions and where N = 2N' 
is the number of spins. In Ref. [H this is called the 
one-dimensional compass model, whereas the authors of 
Ref. [13 reserved that name for the special case with 
Ji = and J2 ~ Li. It can be seen as a ID anisotropic 
XY model with alternating interactions • In Ref. [l^ 
the model in ([T]) was studied with Ji, J2 and Li con- 
strained by a particular dependence on a single param- 
eter a € [0,1], with the construction mirrored in the 
interval a € [1,2]. The solution based on this specific 
choice of interaction parameters indicated a first-order 
QPT at a = 1. However, it was not clear whether this 
transition is intrinsic to the system, or an artifact of the 
singular parameterization of the interactions. To resolve 
this issue, the authors of Ref. [H revisited the problem 
and considered a more general model [isl ]. with the sin- 
gular parameterization removed, 

TV' 

^ = ^\- Jl^2i-1^2i + J2<72i-1^2i 
1=1 

+ Ll(T2i(^2i+l + L2(^2i'^2i+1 ] ; (2) 

with, as before, periodic boundary conditions and N = 
2N' the number of spins. The model was diagonalized ex- 
actly by a direct Jordan- Wigner transformation, yielding 
an integral expression for the ground state energy in the 
thermodynamic limit. It was found that there is a first- 
order QPT when the curves (Ji,L2)(a) or (J2,Li)(q;), 
parameterized by a, p ass through (0, 0), thus confirming 
the finding in Ref. [l2. This is due to a cusp in the energy 
surfaces at (Ji,L2) = (0,0) and (J2,ii) = (0,0), respec- 
tively. In fact, by carefully analyzing the energy surfaces 
obtained in Ref. we find that there is a second-order 
QPT when crossing the lines Ji = L2 or J2 ~ Li, with 
a vanishing excitation energy gap. As these lines pass 
through the point (0, 0) of the first-order QPT, this ex- 
plains the remnants of a second-order QPT that were 
glimpsed in Ref. 

To obtain a transparent picture of the multicriticality 
it is actually more instructive to focus on the ID quantum 
compass model in Eq. ^ [identical to the Hamiltonian 
in when L2 = 0]. Moreover, in order to derive var- 
ious spin correlation functions it is advantageous to use 
the approach developed in Ref. [l^. We here review this 



method, adapted to our case with Ji and J2/L1 in Eq. ([T|) 
being two independent parameters. 

We work in the conventional basis { | T ) , I i ) } of eigen- 
states of the af operators. The only terms in the Hamil- 
tonian ^ that flip spins are those which contain cr^ 
operators, and these terms all have the form 'y2i~i'^2i- 
Thus, the only transitions that can occur are simulta- 
neous flips of spins linked by {2i— l,2i} bonds ("odd 
bonds"). This means that the Hilbert space of the sys- 
tem can be divided into subspaces which are not mixed by 
the Hamiltonian. Each subspace can be labeled by a vec- 
tor s = (si, S2, . . . , s^r/), with element s; = 1 when two 
spins linked by an odd bond {2i — 1, 2i\ are parallel, and 
Si — when they are antiparallel. The Hamiltonian can 
then be diagonalized in each subspace s independently. 
For any subspace s the terms involving crli.iclj will only 
yield a constant contribution Cs(Ji), given by 

N' 

CsiJi) = Ji ^ a^,_,a^, = -MN' - 2s), (3) 

i = i 

n' 

where s = ^h*^ number of parallel odd bonds 

in the subspace s. The Hamiltonian ([1]) then becomes 

N' 

Hs=L,Y,[f ^-i'^^^ - - 2s) . (4) 

Next, define a set of operators 

rf ^ -(in)(iTi + iiT>(ni), (5) 
r'i ^ -(-i)^-^^^(in)(ni-iiT)(iTi), (6) 

for the antiparallel (s^ = 0) odd bond {2i — 1, 2i}, and 

rf ^ -(|TT)(UI + lii)(TTI), (7) 
rf ^ -(-i)^-^^^(ITT>(TT|-|U)(UI), (8) 

for the parallel (s; = 1) odd bond {2i — 1, 2i}. This maps 
the Hamiltonian ([1]) onto 

N'-l 

i=i ^ 

L,[^T^^, + (-l)^r^,Tf ] + C.(Ji) , (9) 

which is the exactly solved p^.[l7j one-dimensional quan- 
tum Ising model with periodic (antiperiodic) boundary 
conditions when s is even (odd) The Hamiltonian 

after a Jordan- Wigner transformation to fermionic c,- op- 
erators and a Fourier transformation to Cj, operators (see 
Ref. m for details), will then be 

= Li^[2((J2/Li)- cos fc)ctc, 
fc 

+ *sinA:(cLfecj; +c_j^Cfe)] 

- J2N' + C'siJi) , (10) 
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where ± denote the separate subspaces of even or odd 
number of c fcrmions, and k takes the values k = 
0, ±j^, ±2j^, . . . , ^ or fc - ±if^, ±1^, . . . , ±i(iV'- 

■'^hen s + X]^i(c]c^) is odd or even, respectively. 

By a Bogoliubov transformation to fermionic 7^, quasi- 
particles, the Hamiltonian in (fTO|) can be expressed on 
diagonal form, 

Hf^J2''^(^hk-l) + CsiJi), (11) 

k 

where efe( J2/L1) = 2Li (1 + ( Js/Li)^ - 2( J2/L1) cos A:)i/2 
and Cs(Ji) = — Ji(7V' — 2s). For s even and a + sub- 
space (or a — subspace with J2/L1 < 1) the number of 
7j, particles is even in that subspace. But for s even and 

a — subspace with J2/L1 > 1, the number of 7^. particles 
is odd. Similarly, for s odd and a — subspace (or a + 
subspace with J2/L1 < 1), there must be an odd num- 
ber of 7j. particles in that subspace. Finally, s odd, + 
subspace, and J2/L1 > 1, means an even number of 7^, 
particles. 

When Ji < 0, the term Cs(Ji) in Eq. ((TT]) forces the 
ground state to be in the s=N' subspace, and not in the 
■3=0 subspace as for Ji > 0. The ground state energy in 
the thermodynamic limit is then given by 

Eo {Ji,J2/Li) = -JiN' - — / dk efc( J2/L1) (12) 

when Ji > 0, and 

Eo{Ji,J2/Li) = JiN'-—J dk€k{J2/Li) (13) 

when Ji < 0. Without loss of generality, we shall restrict 
ourselves to positive values of J2 and ii in what follows. 

III. QUANTUM PHASES AND PHASE 
TRANSITIONS 

The expressions for the ground state energy in 
Eqs. dni) and ^ allow us to find the QPTs of the 
model (U)). There is a discontinuity in the derivative 
dEo/dJi, and therefore a first-order OPT, when passing 
through Ji = 0, as observed in Ref. But we also note 
that there is a second-order QPT when passing through 
J2 = Li, since the second derivatives with respect to 
J2/L1 of the integrals in Eqs. ([T2|) and (fT3|) arc divergent 
at J2/ Li = 1, without any singularity in first derivatives. 
We can thus plot the phase diagram, see Fig. [U which 
shows how the lines of first- and second-order QPTs meet 
at the multicritical point at Ji = 0, J2/L1 = 1. We now 
turn to investigate the different phases separated by these 
transitions. 

Let Ji 7^ 0. From Eq. (HH]) it follows that in the 
ground state (c^^c^) = and (cjcg) = 1 ((cJcq) = 0) 
when J2/L1 < 1 { J2/L1 > 1). We then find that in 
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FIG. 1: (Color online.) Phase diagram of the one-dimensional 
quantum compass model [Eq. ^] . There is a first-order QPT 
at Ji = which separates the phases with (crfi-icli) — 1 
(I. and III.) from those with {<72i-i'^2i) = ~1 (H- and IV.). 
There is also a second-order QPT at J2 /Li = 1 separating the 
phases with (erf) = (I. and II.) from those with (af ) 7^ 
(III. and IV.). The dashed lines show the three paths (a), (b) 
and (c) that are used in Fig. (2] All paths start at the point 
Ji/Li = 1, J2/Z/1 = 0, where the Hamiltonian ([1} reduces to 
that of the one-dimensional quantum Ising model. 



the thermodynamic limit, for J2/L1 < 1 the two states 
with no Bogoliubov quasiparticlcs, s = 0, and either in 
the subspace -I- or — will be ground states, i.e. double 
degeneracy. For J2/L1 > 1, only the state in the + sub- 
space will be a ground state, i.e. no degeneracy. The 
degeneracy for J2/L1 < 1 will be lifted in a finite sys- 
tem since then the summation over k in Eq. (|11[) will be 
different in the ± subspaces. 

When Ji = the energy eigenvalues are independent of 
s, since (7^(0) = 0. Thus for a finite system, when Ji = 
and J2/L1 7^ 1, the ground state degeneracy d will hcd = 
2N'-1 ^ 2^/2-1 , since each state with no jk quasiparticle 
in every {s, -I-} subspace with even s will have the lowest 
energy. When Ji = and J2/L1 = 1, the 7g quasiparticle 
gets gapless, so that the state with one 70 quasiparticle 
in every {s, +} subspace with odd s also will have the 
lowest energy. Then d becomes d = 2^ = 2^/^ for a 
finite system. In the thermodynamic limit each subspace 
s with s even contributes two states to the ground state 
degeneracy at Ji = 0, J2/L1 = 1 (the states without 7fc 
quasiparticles in both ± subspaces), and each subspace s 
with s odd also contributes two states [the state with one 
7q quasiparticle in the + subspace and the state with one 
^n/N' (which becomes 7o when N' —^ 00) quasiparticle in 
the — subspace] . Therefore the ground state degeneracy d 
will be 2 X 2^/^ when N ^ 00. We thus confirm the result 
for the degeneracy at the transition point a = 1 (i.e. Ji = 
and J2/L1 = 1) derived in Ref. [l^. The discrepancy 
with the result in Ref. [isl where the degeneracy was found 
to be d = 2^/2-1 remains unexplained. 



4 



(a) J^/L^ =2(1 -J^/L^] 
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(b) J2/L^ = 1 - J/L^ 
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(c) J2/L^ = (1/2)(1 - J^/L^ 
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FIG. 2: (Color online.) Summary of the properties of the one-dimensional compass model [Eq. ([TJ] in the thermodynamic limit. 
The plots are along the three paths (a) J2/L1 = 2(1 - Ji/Li), (b) J2/L1 = (1 - Ji/Li) and (c) J2/L1 = (1/2)(1 - Ji/Li), 
shown in the phase diagram in Fig. [l] Left column: path (a). Middle column: path (b). Right column: path (c). Upper row: 
energy gap A. Middle row: two-point functions in the ground state. (Blue line (1): (crfi-icrfj), red line (2): (cr|iCrfj+i), blue 
line with bars (3): (o"2iO"2i+i), red line with circles (4): (crli-icfi))- Lower row: ground state degeneracy d. 



The energy gap between the ground state and the low- 
est excitations is obtained from the diagonalized Hamiho- 
nian in Eq. pip and the subsequent rules for the number 
of 7j. particles allowed in each subspace. In the ther- 
modynamic limit, the gaps are Ai = 4Li|l — J2/Li\ for 
the state in the s = subspace with one 7q and one 
^±2tt/n' particle, A2 = 2|Ji| for the state in the s = 1 
subspace with no 7^, particles (this is only possible when 
^2/^1 > 1), A3 = 4|Ji| for the state in the s = 2 sub- 
space with no 7^ particles, and A4 = 2Li|l — J2/Li\ 

for the state in the 5 = subspace with one 7q particle 
(which is only possible when J2/L1 > 1). At each point 
in the phase diagram the excitation energy gap from the 
ground state is now given by the smallest of Ai^2,3,4- This 
gap is plotted in Fig. [51 along three representative lines 
in the phase diagram (cf . Fig. [T]) . The first- and second- 
order QPTs are marked by vanishing energy gaps. 

As previously seen, the ground state is in the s=N' 
subspace when Ji < 0, and in the s—0 subspace when 
Ji > 0. Therefore the ground state two-point functions 
(^li-i'^'li) are discontinuous at Ji = 0, with (cr|i_icr|i) = 
1 when Ji < and (cfi-i^'li) = ~1 when Ji > 0. On the 
other hand, {o'2i'^2i+i) = in all phases, since the oper- 
ator a-2i'^2i+i takes the ground state out of its subspace 
s. The mapping to the QIM [Eqs. (O to ^] imphes 



that cr|i_icrf, ^ -rf and cr|,cr|i+i ^ -t^t^+i, so that 
K.-iC^I,) = -(^f > and = -(rfrf+i). 

The characteristics of the different phases and phase 
transitions of the model are graphically summarized in 
Fig. [21 where we plot the energy gap, the two-point func- 
tions and the ground state degeneracy along the three 
paths (a), (b) and (c) shown in Fig. [2 It is seen in Fig. [21 
that the continuous QPT at J2IL1 = 1 is associated with 
a continuous transition from dominating antiparallel or- 
dering of spin z components on even bonds for J2/ Li < 1, 
to dominating antiparallel ordering of spin x components 
on odd bonds for J2/L1 > 1. The order parameter (tTji) 
or {o'2i~i) only non-zero in the ordered phase given by 
J2/L1 < 1, since = \{ai,_,)\ = (rf). The first- 

order QPT at Ji = corresponds to a = 1 in Ref. [l^- 
We see that the mixed first- and second-order features 
of the QPT in the one-parameter model studied by the 
authors of Ref. [l^ come about because of their param- 
eterization, where the multicritical point is approached 
along a path [(b) in Fig. [T] with projections along both 
of the transition lines. 

At Ji = the Hamiltonian (|T]) reduces to what Ref. [T^ 
refers to as the proper one-dimensional QCM, 

N' 

= E [ •^2^2,-1^ + ] • (14) 

i=l 



5 



It was shown above that the ground state degeneracy for 
a finite system is d = 2^/^~^, except at the second-order 
QPT at J2 ~ Li, where it is d = 2^/^. The large ground 
state degeneracy can thus be easily understood as a con- 
sequence of the model ([M]) being at the level crossings of 
the first-order QPT of the extended ID QCM (P). 

IV. ENTANGLEMENT 

Given the exact solution of the ID QCM, we have a 
rare opportunity to analytically probe for the entangle- 
ment in the ground state of a highly complex system of 
coupled qubits. We here focus on two of the most fre- 
quently used entanglement measures: concurrence and 
block entanglement entropy |19|] . 

A. Concurrence 

The concurrence of two spins at sites i and j is ob- 
tained from their reduced density matrix pij, which can 
be expanded as [l^] 

P^^ = lY.'^»t''s- (15) 

where — Q,x,y,z, with = 1^. Note that the 
ground state expectation values of af , erf, erf cr| and ctJ 
(and i <-> j) are zero, since these operators flip only one 
spin on an odd bond, giving a state in a different subspace 
than the ground state. The correlation function (erf ctJ) 
must also be zero, since the matrix erf ctJ is imaginary and 
Pij must be real as the Hamiltonian is real. Therefore, 
the reduced density matrix p5p becomes 

+ {<r\ay)ayay + {ala^^)ata^,). (16) 

The concurrence C{pij) is now given by [2l| C{pij) = 
max{ , Ai — A2 — A3 — A4 }, where Ai > A2 > A3 > A4 
are the non-negative real eigenvalues of the Hcrmitian 
matrix R = Pij Pij Pij ■ In this expression, = 
{afaj)p*j {crfaj), where p*j is the complex conjugate of 
Pij in the given basis. Here, the eigenvalues ^'1^2, 3, 4 of R 
are (without ordering) 




xv'l-«) + (^j)-«^|> 
±\{afa]) + {ayay)\\. (18) 



It is immediately clear that the concurrence of two spins 
that are not on the same odd bond is zero, since then 
(erferj) = (er^'erj) ~ which gives ri ~ r2 and r^ = r4. 
For two spins that are on the same odd bond {2i — 1, 2i}, 
we have that (erf^.ierfi) = —{rf)- When Ji > 0, the 
ground state is in the subspace 5=0, so then (er|j_]^) = 
-(f^l»> = "(t/) and (er|,_iCT|,) = -1. We also have 
that (cr|,_ier|,) = (o-2j-iC^2»>: since the matrix al^_^al^ 
is the same as erfj_]^er|j in the subspace s=0. Then 
g(p2»-i 2») = ma x{0 , r3 - r^} = -(erf,;_icrf,) = (rf), 
since yl — (rfP > {Tf)- On the other hand, when 
Ji < 0, the ground state is in the subspace s=N', 
where (a^.^af,) - -(rf), (a|,_i) = (a|,) = -(rf), 
and (er|j_ier2j) = 1. Also, in the subspace s=N' the 
matrix er|j_j^er|- is the same as — crfj.j^erjj, so that now 
(er|._^er|-) = — {o'2i-i'^2i) ■ ^6 then get the result that 
C(p2i-i,2») = max{0, ri - rj} = -(erf^.^erf,) (rf). 
Thus, for all values of Ji, the only pairs of spins with 
non-zero concurrence are those on the same odd bond, 
and for them C{p2i-i.2i) ~ (^f). The pairwise concur- 
rence therefore only depends on the parameter J2/L1, 
which plays the role of the magnetic field in the QIM 
From the well-known solution [H, [l3| of the ID QIM, it 
follows that in the thermodynamic limit the concurrence 
has a diverging first derivative across the curve J2/L1 = 1 
of the second-order QPT. Note that the first-order QPT 
at Ji = is not signaled by the pairwise concurrence. In 
particular, this means that when going through the mul- 
ticritical point, the concurrence will behave as if it was 
a pure second-order transition, despite the transition ac- 
tually being of first order. 



B. Block entropy 

The block entropy is the von Neumann entropy of a 
subsystem consisting of an entire block of adjacent spins, 
thus giving a measure of the amount of entanglement 
between the block and the rest of the system [1^. We 
will now show that a block of an even number of er spins 
in the ID QCM ([T]) that fully covers an integer number 
of odd bonds will have the same block entanglement as 
half the number of r spins in the corresponding QIM ([9]). 
Let us first consider the case when Ji > 0. 

The reduced density matrix p^ for a block of L spins 
eri, . . . , er^, where L is even and eri and er2 arc on the 
same odd bond, can be expanded as [l^l 

P^ = ^ H (19) 

fJ.l,---,fJ.L 

where /^i, . . . , /i^ are summed over 0, x, y, z. For the ex- 
pectation value (er^'^ ...cr^^) not to be zero, the only 
allowed pairs of er operators on every odd bond in the 



6 



summation are 
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Thus the reduced density matrix (fT9|) becomes 

l'l,...,I/l,/2 

where vi, ... ,1^1^/2 ^re summed over 0,x,y,z. This is 
precisely the reduced density matrix of L/2 spins in the 
QIM. The block entropy Sl is then the same as the block 
entropy of L/2 spins in the QIM with effective trans- 
verse field equal to J2/L1. This will apply equally well 
also when Ji < 0. The phase transition at J2/L1 ~ 1 
is therefore in the Ising universality class, since this is 
uniquely determined by the scaling of the block entropy 
at criticality The first-order QPT at Ji = and 

with J2/L1 ^ 1, corresponds to a non-critical QIM, with 
saturated block entropy Sl when L —f 00. Its value will 
be the same whether Ji 0^ or Ji ^ 0+. 

V. CONCLUSIONS 

To summarize, we have performed an exact analytical 
study of the ID quantum compass model, using a map- 
ping to the quantum Ising model, following the approach 
in Ref. We identify four distinct ground state phases. 



separated by two intersecting transition lines. One of 
these defines a line of second-order Ising-like transitions, 
while the other is a line of first-order transitions (cf. 
Fig. [T]). The point of intersection, where the first-order 
quantum phase transition identified by Brzezicki et al. 
[12| takes place, thus defines a multicritical point. This 
explains the apparently exotic behavior at the transition 
found by these authors. In particular, the appearance 
of a diverging correlation length for certain spin corre- 
lations finds a natural explanation once the multicriti- 
cality of the transition point has been recognized. One 
may ask whether the quantum phase transition in the 2D 
quantum compass model maybe also plays out at a mul- 
ticritical point, in analogy to the ID model? If so, this 
could possibly explain the notorious difficulty in identi- 
fying the character of the transition, as evidenced by the 
conflicting results in Ref. 9 and Refs. d, [13, fill . 

Our results for the entanglement show that the only 
effect on the ground state when going through the first- 
order transitions is that a correlation function for neigh- 
boring spins on odd bonds changes sign, without any ef- 
fect on the entanglement measures studied. First-order 
QPTs are generally associated with a discontinuity in 
concurrence, but "accidental" exceptions to this rule are 
possible [2^. We here have an example thereof. 
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